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Abstract 

Recently, M.K.-H. Kiessling and A.S. Tahvildar-Zadeh proved that 
a unique global classical solution to the relativistic Vlasov-Poisson sys- 
tem exists whenever the positive, integrable initial datum is spheri- 
cally symmetric, compactly supported in momentum space, vanishes 
on characteristics with vanishing angular momentum, and for j3 > 3/2 
has & -norm strictly below a positive, critical value Cp. Everything 
else being equal, data leading to finite time blow-up can be found 
with £^-norm surpassing Cp for any /3 > 1, with Cp > if and only 
if f3 > 3/2. In their paper, the critical value for (3 = 3/2 is calcu- 
lated explicitly while the value for all other f3 is merely characterized 
as the infimum of a functional over an appropriate function space. 
In this work, the existence of minimizers is established, and the ex- 
act expression of Cp is calculated in terms of the famous Lane-Emden 
functions. Numerical computations of the Cp are presented along with 
some elementary asymptotics near the critical exponent 3/2. 
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1 Introduction 



The relativistic Vlasov-Poisson (rVP) system is given by 



d t + 



p 



= -V q ±V q <p t (q) 



V p / t (p,?) = 



rVP 



± . 



< 



4nff t (p,q) d 3 p 



<Pt(q) x -\q 



i 



as \q\ — > oo; 



v 



rVP + models a system with repulsive interaction while rVP~ models a sys- 
tem with attractive interaction. One of the earliest papers to appear on the 
subject is |GS85] wherein Glassey and Schaeffer show that global classical 
solutions to rVP ± will exist for initial data that are spherically symmetric, 
compactly supported in momentum space, vanish on characteristics with van- 
ishing angular momentum, and have £°°-norm below a critical constant 
with C+ = oo and < oo. More recently, Hadzic and Rein ( [H R07] ) showed 
the non-linear stability of a wide class of steady-state solutions of rVP - 
against certain allowable perturbations utilizing energy- Casimir functionals. 
Shortly thereafter, Lemou, Mehats, and Raphael f |LMR08al ILMR09| ) in- 
vestigated non-linear stability and the formation of singularities in rVP - 
through concentration compactness techniques. 

In this work, we focus exclusively on the attractive case and henceforth 
suppress the superscript on both rVP~ and C^. In |KTZ 08j. Kiessling and 
Tahvildar-Zadeh prove that a unique global classical solution to the rela- 
tivistic Vlasov-Poisson system exists whenever the positive, integrable initial 
datum /o is spherically symmetric, compactly supported in momentum space, 
vanishes on characteristics with vanishing angular momentum, and has £ /3 - 
norm below a critical constant Cp, with Cp > if and only if /3 > 3/2. The 
constant Cp is critical in the sense that, everything else being equal, initial 
data can be found which lead to blow-up in finite time if their £^-norm is 
allowed to be ever so slightly bigger than Cp. This critical constant is given 
by the following minimization problem: 




(1) 
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where D £^ denotes the set of probability measures on M 6 with finite first 
moment which are absolutely continuous with respect to Lebesgue measure 
having density in £ /3 . The functionals and S q are given by 



S;(f) = JJ\p\f( Pl q)d 3 pd 3 q 




\q-q'\ 



d 6 p' d 6 p d 6 q' d 6 q. 



(3) 
(4) 



It is further shown that 



3V 15 
8/ 16 



<C B < 



45 



87T5 

nLi(^ + 2/3)r(/3 + l) 



(5) 



In particular, the optimal constant for (3 = | is explicitly calculated in the 

paper. Its value is | (jf) 3 , which is equal to the one given by either bound. 

In this paper, we address the determination of the critical norm for all 
remaining cases. We begin in the next section by proving the existence of 
minimizers for $^ in a slightly larger class of functions Qp to be defined 
below. After that, we characterize the minimizers via variational techniques 
and find that they are the well known Lane-Emden polytropes. This will 
show that the minimizers are actually in our original space n ff. Finally, 
we compute the optimal constant in terms of the parameter (3 and the first 
zero and corresponding slope of the standard polytropes. At the end, we 
mention some numerical results pertaining to the calculation of Cp. 



2 Existence of Minimizers 

To begin with, as remarked in the "note added" in |KTZ08j . the variational 
problem (1),(2) is equivalent after rescaling to the one given by Lemou, 
Mehats, and Raphael in |LMR08b] . who designed it for different purposes, 
namely to study blow up dynamics for the relativistic Vlasov-Poisson sys- 
tem. Instead of giving a detailed analysis, they refer the reader to their 
earlier work |LMR08a] where they study an analogous variational principle 
(given in formula 1.18) for the non-relativistic Vlasov-Poisson system. Thus, 
in principle we could build on their results to evaluate Cp. Instead, here we 
first give a somewhat different proof of the existence and characterization of 
minimizers by combining the techniques of Weinstein ([W83]) (also referred 
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to by Lemou, Mehats, and Raphael) with those of Lieb and Simon f [LS77j ). 
This strategy has certain advantages, as we shall see in the next section. 

We find it convenient (for reasons explained below) to expand the class 
of functions over which we attempt to minimize <3>^. To that end, define 

Sip = {/ : R 6 -> R : / > 0, \\f\U + ||b[/||i + \\f\\ p < 00} 

and note that ^1 fl £^ C Qp (indeed, a function / e ftp will also be in 
^Pi PI £^ whenever |||?|/||i is finite and the £ 1 -norm is equal to 1). Since we 
have allowed functions of arbitrary £ 1 -norm into our considerations, we need 
to adjust our definition of <&g. Assuming > and inserting y^U into 
$3, we arrive at the appropriate functional: 

/ S u (f) \ 3(1 ~^ 3 _i 

Mf) ^\^U)J 11/11/3 11/111 (6) 

We now seek to minimize Qp over this enlarged space of functions. In the 
next section, we will show that there exist minimizers for $g over Qp that 
arc also in ^n^. 

The demonstration that minimizers exist closely follows Weinstein ( |W83j ). 
Since $>p{f) > for all functions in Qp, we can find a minimizing sequence 
{fp, n } C Qp so that 

Cp ee inf = lim $p(fp, n ). (7) 

Since we are minimizing over a larger class of functions, Cp < Cp and so the 
upper bound noted above still holds. Applying the well-known procedure of 
spherically symmetric equi-measurable rearrangements (c.f. |LL01l Chapter 
3]), we can assume fp, n (p,q) = fp,n(\p\, \q\,&) (committing a slight abuse of 
notation) where 9 is the angle between p and q. 

For any positive real numbers k, A and fi, the triple family of scaling 

f*,\,i*(p,q) = V /(Ap,«g) (8) 

leaves $^ invariant while 

£p (//c,A,/i) : 

^ij(/k,A,/j) : 
H/rejAj/i || r 





(9) 




(10) 


^ 3 11/11,- 


(11) 
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Taking advantage of this scaling invariance, we can assume that our mini- 
mizing sequence has the following properties: 



fp, n (p,q) = /An(bl,M>0), (12) 

£;{U,n) = 1, (13) 
Wfanh = 1, (14) 
Wfonb = I" (15) 

Hence, the Banach-Alaoglu theorem and the reflexivity of £^ for j3 G (1, oo) 
give us a function fp G £^ such that some subsequence of {fp, n }nLi converges 
weakly in £^ to Z^. Without loss of generality we assume that we have 
already extracted this subsequence. Standard arguments concerning lower 
semi-continuity and weak convergence (c.f. |FL07| ) show that 

W/0 < 1> (16) 

< i, (17) 

ll//»ll/J < 1, (18) 

so that fp G flp (this is the advantage of expanding the space of functions as 
showing that fp G fl £^ is difficult at this point). Thus, we can conclude 
the following is true: 

C P <hifp) =(5^)) 3(1 "" ) Wfeh MpC* <(-^(//»))" s(1 "* ) . (19) 

Since by construction 

hm ^(J^) = lim (_f (/^ ))- 3 d-?) = fy, (20) 

we will be done if we can show that £ q (fp tn ) converges to £ q (fp). Unfortu- 
nately, S q is upper semi-continuous with respect to weak convergence, and 
so all we can immediately conclude is that 

£ q (fp) > lim E q (fp in ). (21) 

To show the convergence of £ q (fp tn ) to £ q (fp), we first rewrite the poten- 
tial energy functional as 



W) = -\j Pfil) K f {q) d 3 q, (22) 
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where 



and 



p f (q) = J f(p,q)d 3 P , (23) 

K f (q) = Pf * lldl" 1 ^) = / d\>. (24) 

J \q — q'\ 

We seek bounds on pf which will then imply bounds on Kf ffn . To that 
end, Lemma 4.3 in [KTZ08J implies that 

Wpfh^cmfV^if) 1 -" (25) 

with exponents given by 

1b = —t. and 713 — — . (2d) 

w 3(3-2 ^ 4/3-3 K ' 

Note that 7^ is an increasing function of (3 and that the limiting case of 
(3 = § gives 73 = f. Thus p fp and p fp<n are in ^(M 3 ) n £ 7 ^(R 3 ) for all 
(3 and all n. Also note that the sequence {pfg n }^=i is uniformly bounded 
in £ 7 ^-norm, and so some subsequence must converge weakly in this space. 
Standard arguments show that this weak limit must equal pf p a.e. 

Next, we can conclude that Kf and Kf are in £^ C (M 3 ) for 3 < a < 
f fLLOil Theorem 10.2]). Note that ±^ is also an increasing function 
of (3 and the limiting case of (3 = | makes this exponent equal to 6. We 
can turn the local estimates into global ones for all a > 3 via the following 
growth estimate on the potential. 

For any spherically symmetric f(p,q), the marginal mass distribution 
Pf{q) will also be spherically symmetric. Hence, the well known formula for 
the potential of a spherically symmetric mass distribution gives 

47T f\ q \ f°° 

K f (\q\) = y-r / p f (r) r 2 dr + An p f (r)rdr (27) 

\q\ JO J\q\ 

Air Air f°° 

< — - / pf(r) r 2 dr + — pf{r) r 2 dr (28) 

\q\ Jo \q\ J\q\ 

- w (29) 

where in the last step we have used that for our purposes ||p/||i < 1. So, if 
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K f e ££ C (M 3 ) for a > 3, then we have for any R > 

[KJ(q)d 3 q = [ Kf(q)d 3 q+f Kf(q) d 3 q (30) 
J J\q\<R J\q\>R 

£ L K ^ d3 « + (^w^- (3i) 

Thus, we may conclude that for the distributions we are considering, Kf„ n 
and K ffj are in £ a (M 3 ) for 3 < a < m ^ ± . Note that this proves that p f K f 
is indeed in £ 1 (R 3 ) for spherically symmetric / G Qp. 

As with Pf fjn , we have a bound on the £ a -norm of Kf„ n that is indepen- 
dent of n, and so some subsequence of {^/«„}^Li must converge weakly in 
this space. As usual, this weak limit must equal Kf p a.e. 

We can actually say something much stronger about the convergence of 
{Kj p n } to Kf - namely that any subsequence converging weakly to Kf p will 
actually converge strongly in £ r on sets of finite measure for any r < 12 . 
To see this, we note that since V • VUffg n is in 2} fl £ 7 ' 3 , we know that 
VK ff) n is locally in fit n £ K ^ where «/j = . Thus, Vify, _„ is locally in 
£ 2 whenever /3 > |. Arguments like the one leading to fl29|) for if/g n show 
that these functions decay rapidly enough at infinity to be in £ 2 proper, and 
that the decay is independent of n. Again, some subsequence of {VKf pn \ 
will converge weakly in £ 2 and must converge to VKf . Following the proof 
of [LLOlj. Theorem 8.6], we have the strong convergence stated above. 

Combining our upper bound for Kf and the local strong convergence gives 
us that {Kf n } converges strongly to Kf p in £ a for 3 < a < 12/3 ~ 9 : 

\\Kf,. n -K f X = \\{K f ,. n -K u )xB R( o)\\a+ [ \K f ^-K ff rd 3 q(32) 

Jb r (o) 

< \\{K U , n -K f0 )xB R{O )\\ a a + [ j^-d 3 q (33) 

Jbuo) \1\ 



< II (K fp<n - K fp ) X B R (o)\\ a a + Q , „ a _ 3 - (34) 



'B%(0) |?| 
47T 

(a-~3)i2 



Finally, we are in a position to show the convergence of 8 q (fp in ) to 8 q (fp). 
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We assume that > | and hence pf G £s and ii'/ G £ 6 . 



2\£q(fp)—£q{fj3,r. 



< 



Pf?,M(K u Jq)-K ffi (q))d 3 q 



(35) 



+ 



(pfpM - PfM)) K fM) d *<i 



< Wpu 



\ K fp,n- K fph 



+ 



(36) 



(37) 



The first term in the last inequality can be made arbitrarily small since Kf„ n 

converges strongly to Kf p in £ 6 for f3 > | and the £§-norm of pf p n is bounded 
above independently of n by standard interpolation estimates. The second 
term can be made arbitrarily small by the weak convergence of p/ (3n to pf 
in £s (after possibly another subsequence extraction) and the fact that Kf p 
is in the dual space - £ 6 . 

We see that inequality (fl9|) is saturated, and we have that 



Note that this forces 



\\M\i 

1/3 



1, 
1, 
1. 



(38) 



(39) 
(40) 
(41) 



Thus, we have shown the existence of minimizers for $^ over the space Qp 
satisfying the properties given above. 

In the next section, we show that fp can be chosen so that it is compactly 
supported for (3 > |, so that among all possible minimizers for $^ there is 
indeed one in D £^ - proving that Cp = Cp. 



3 Identification of Minimizers 



We want to find the infimum of §p over the space of functions flp introduced 
in the last section. Following an idea in Lieb-Simon ( |LS77j ) we first note 
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that £lp is a convex space of functions so that, if fp is a minimizer of our 
functional over Qp and rj is any function in this space, then for all < t < 1 
we have that (1 — t)fp + ti] e f|g. Consequently, we can consider 



t=o+ 



^((l-t^ + ty), 



(42) 



where by t = + we have in mind the one-sided Gateaux derivative from the 
right at zero. This technique avoids the difficulty that arbitrary variations 
of a given / e ftp may become negative (and hence no longer belong to Qp). 
Direct calculation gives us that this derivative is 

II//3II/3 ^ ll/plll J 

We now compute the indicated derivatives separately. We begin with the 
ultra-relativistic kinetic energy: 



d 
dt 



£;((l-t)fp + trj) 

t=o+ 

Next, we easily compute: 
d 



-£p(fp) + J J \p\ V(P, q)d 3 pd 3 q. (44) 



dt 



t)U + tv\\r = 



t=0+ 



WfiiW 



1-r 



-\\M\ T r + J J (MP, 



q)) r 1 r](p,q) d 3 pd 3 q 



(45) 



Finally, we find: 
d 



j t J q ((l- t )fp + trj) =-2(S q (fp))-JjKp(q) V (p,q)d 3 pd 3 q(A6) 

where we have used Kp as in the previous section. Inserting these into our 
derivative above, collecting terms and noting that all the constant terms 
cancel (i.e. those terms not involving an integration against rf) yields 



d_ 

dt 



^((l-t)fp + tr]) 



t=o+ 



3- 



\p\ Kp{q) \ | (frM)'- 1 | 2-| 



(47) 

r)(p, q) d z p d 3 q. 



9 



Since the indicator function of any set of finite measure is in Qp, we are 
tempted to simply conclude that 



3- 



\P\ 



Kp{q) \ , (fp(p,q)) 



or equivalently that 



2 - 1 



= o 



(4£ 



//J(P,9) = II//J||J 



0-1 



P 



\P\ 



KM 

-Wp) %<Jfi)J \\M\ 



2 - ^ 



p-i 



(49) 



Such a function cannot be in Qp since whenever 



3 - 



Kp{q) 



\P\ 



(3 J \-s q (fp) s;(fp) 



2-4 



< o, 



(50) 



we get complex values for fp in general. Hence, we take 

Kp(q) \p\ \ 2- 



//3(P,?) = 



g 
0-1 

P 



3-1 



,(51) 



where (■)+ means the positive part of the argument. 

Since we have altered the natural minimizer so that it lies in Qp, we need 
to examine the effect on the one-sided Gateaux derivatives of our functional. 
Let A be the support of our minimizer fp. Every r\ G ftp can be decomposed 
as 77 = t]Xa + VXa c - Hence, 



d 
dt 



$p((l-t)fp + tr,) 



t=o+ 



A? 



3-1 



\p\ k p(q) \ , (fp(p,q)) 



Pi \^(fp) -Wp) 



r-i + 2-| 



\\fp\\ 



(52) 

rj(p, q) d 3 p d 3 q. 



But on the the set A c , we have that fp = and that the integrand is strictly 
positive. Hence, for any r] 6 ftp 



d_ 

dt 



$p{(l-t)fp + t V )>0, 



(53) 



t=o+ 



showing that fp as defined in floTT) is indeed a minimizer for $p over 
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We now use scaling invariance to pick out a special minimizer that we 
will use to calculate Cp. To that end, recall that the triple family of scalings 
as given in (jSj) 

/«,a,/*(p>?) = H f(\p, Kq) 

(where k, A and /i are positive real numbers) leaves Qp invariant. As noted in 
the previous section, using this scaling we can choose fp so that its £ 1 -norm is 
1. We break from the remaining choices above by requiring —£ q (fp) = £p(fp) 
(as opposed to £p(fp) = 1). Hence, we can no longer choose the £^-norm to 
be 1. Such choices are completely arbitrary since we have an infinite number 
of minimizers from which to pick. However, these choices give a minimizer 
of particularly nice form. 
Our choices so far force 

^=\\M\l-^jf-y (54) 



and 



i 

which leaves us some room to choose either k or A in a convenient way. Using 
(IBTj) . we see that 



K = ^ir = \\M\ 1 [^ m ) , (55) 



fp-, K ,x,n(p,q) = H 



3V, n pfKp( Kq ) \Xp\ \ {'~-p 



0-1 



(56) 

2_lA ii /ji/3\ 



(57) 



where we have used Ka' x,tJ, (q) = J J — K ' x ^f ' g ^ d 3 p' d 3 q'. If we choose for k 



and A: 



3 V 1 (-We))**-* 



A 1 : *J iimi? w ' (59) 



ii 



then we get the following: 

i 

U,K,\,n(p,q) = {<Pp-,k,\,M) - \p\)+' , ( 60 ) 

where 

ft .,x,M = jj U ^ X q ^ q ' ] d 3 p> d s q' - x, (61) 
and where the constant x can be determined from the constraint 

II/#«,a,mI!i = 1- (62) 

From this point forward, we drop the subscripts k, A, and fi as we now focus 
our attention on this particular minimizer. 

Next, we determine (frp in more detail. We begin by computing the 
marginal mass distribution over configuration space: 

Pf)(q) = J Mp,q) d 3 p, (63) 

= 4tt / ^ q )-\ p \)— \ p \*d\p\, (64) 
Jo 

8tt(/3- l) 3 tl . 

(Ml))t , (65) 



0(20 - 1)(30 - 2) 

where the last line follows by successive integration by parts. By definition, 

Mq) = J J " J fz7f d3 p' d ^ ~ x ( 66 ) 



3;3-2 



- if r (M<t))r 



tmmi (67) 

J \q-q\ 



0(20-l)(30-2) 
so upon taking the negative Laplacian of both sides we have 

3/3-2 

-AMQ) = c{0)(MQ))^ (68) 

with 

= 32^ - If 
lW -/3(2/3-l)(3/)-2)' { ' 

Our previous arguments with spherically symmetric equi-measurable rear- 
rangements show that we can assume (ftp is spherically symmetric (where we 
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take the center of the distribution to be the origin of our coordinate system). 
We write 

MQ) = M\9\) = Mr)- (70) 

Let R/3 be the first zero of 4>p(r). Then on the ball \q\ < Rp the partial 
differential equation above becomes 

(0) = 0, (71) 



Note that the first boundary condition is forced by the differential equation if 
we are to have a finite solution at the origin. The second boundary condition 
will ensure that fp integrates to one (which also specifies the constant x 
which we have conveniently absorbed into the definition of <j)p). 

Noting that the equation fl68|) essentially gives the gravitational potential 
of our mass distribution (up to a sign difference), outside the ball \q\ < Rp 
we must have 

Mq) = i^r - -b- ^ 

\Q\ Rp 

We note that as a potential function, <pp is zero at the boundary of the mass 
distribution and not at infinity (as is usually the case). This fact is easy to 
forget and can be the source of many headaches! 

Since for each j5 > | the associated density fp has compact support (c.f. 
section 5), we can conclude that fp G R & as promised in the previous 
section. Hence, we conclude that Cp (the infimum over fp G is equal 

to the infimum over the larger space Qp, and both are given by Qp(fp). Note 
that this is also true for the critical case (3 = | as shown in [KTZ08J. 



4 Calculation of Cp 

We begin by recalling the calculation above (c.f. f )65|) ) for the the spatial mass 
density associated to fp (that is, the marginal distribution over g-space): 

87r(/3-l) 3 m=2 

p m = m _ 1)m _ 2) (M«))r ■ 
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4.1 A Most Useful Identity 

We begin with a simple integration by parts on the defining PDE for 4>p\ 

c(JS) ! (Ml)) 4 -^ d 3 q = I \Vcf>p{q)\ 2 d\ (73) 

since <f>p(q) = when \q\ = Rp. 

We pair this with the Pohozaev identity: 

[ \VMQ)\ 2 d 3 q = jf\c{P) f (MO))** A - ^. (74) 
JB Rfj (o) 4/} -3 Jb R0 {o) Rp 

This identity can be seen by first noting 

c(/3) / (q ■ VMQ)) (MQ))^ A = [ (9 • ^M<l)) (-A^(g)) d 3 q. 

(75) 

The left-hand integral is fairly easy: 

/ (q-VMl))(M<l))^d 3 q=^^ [ {4>p{q))^d\ (76) 

The right-hand is more involved, and after a lengthy calculation yields: 
/ (q-VMQ))(-^MQ)) A = ~ / |V0,(g)| 2 rf 3 g-^. (77) 

•^(0) Z JB Rp (0) Up 

Finally, combining these two expressions for the Dirichlet integral gives 

/ |Vfc( 9 )|»A = £(^§). P8) 

•/Bi^(O) ^ V^P-3/ 

or equivalent ly 

cW J <*,<»)>F A - | (^|) • (™) 

As we show below, this identity makes computation of the functionals com- 
prising §p remarkably easy! 
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4.2 The & Norm 



rr4>p(q)+ 8 
4tt JJ (M<l) + -\P\)~ \p\ 2 d\p\d 3 q 

87r(/3 - l) 2 /•/•^(«)+ , 3/3-2 



(20-l)(30-2)JJ o 
8tt(/3-1) 3 



(2/3-l)(3/3-2)(4/3-3) 
1 



rr<Pf){q)+ 3g _ 2 
J J (M<l) + -\P\)— d\ P \d 3 q 

(MQ)+)^ d 3 q 



The Ultrarelativistic Kinetic Energy 

\p\ff}(p,q) d 3 p d 3 q 



4vr JJ (Mq) + -\P\)~ \p?d\p\d 3 q 

^^/f (?)+ (^) + -bl)- \P? d\p\ d 3 q 

24tt(/3- l) 3 /yM<?)+ ; 3 ^-2 



m-im-*)lh iM9)+ ~ M) rf|p| 4 

/ 4/3-3 „ 



,8(20 - l)(3/3 - 2)(4/3 - 3) 



4tt(4/3 - 3) 

3/3-3 
^(2/3 - 3) 
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4.4 The Potential Energy 

Our choice of scaling is meant to ensure that the potential energy of our 
minimizers equals their ultra-relativistic kinetic energy. To check this, we 
calculate the potential energy directly. 

- sm = \ //// i^MpA h d s p H d s q (94) 

= - // P-MPM £q (95) 
2 J J \q~ q'\ 

(P - 1) 3 ff (MqU)^ \'Mq') «. « 



m-iw-2)JJ dqdq (96) 



(B-l) 3 f 30-2 



0(20-l)(S0-2) J Br 



\'Mq') d s q , 



d 3 q. 
(97) 



We first work on the bracketed integral before tackling the entire expression: 

/ = ¥w 1 *' = / v a ,« 9 ')-v,(^)dV 

JB Rp (p) 19-91 Jbr p (0) VI?- 9 1/ 

- / (uT^r) V «"W) ■ d °' (98) 
JdB R p(o) \\q-q\J 

^(o) VI9-9I/ 

JdB Rfj (o) \\q-q\/ 

- [ (r-^)v q >Mq')- da' (99) 

JdB R p(o) \\q-q\J 

= 47r^ /J (9) + ^. (100) 
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Inserting this expression into the functional above gives: 



4tt(/3 — 1) £ 



0(2/9-1X3)9-2) 7^ {0) 



3/3-2 

(Mi)) ^ 



Mq)+-b- 



c(/3) 



8vr 

1 Z4/3-3 



(0/3 (g))^ rf 3 g+ ' 



2R p V.2/3-3 

3/3-3 
^(2/3-3)' 



+ 



2i?. 



2i? 



-/b^co) 



d 3 g(101) 

P p{q) d\l02) 

(103) 
(104) 



4.5 The Formula for Cp 

We see from our work above that the ultra-relativistic kinetic energy and 
potential energy are indeed equal for fp (as we chose in our scaling). Thus, 
the only contribution to Cp comes from the £^ norm. Hence, we arrive at 
the following formula: 



C, 







(105) 



5 C/3 and the Standard Lane-Emden 

Polytropes 

Though the formula for Cp given above is rather elegant, Rp is only defined 
implicity by our requirement that \\fp\\i = 1. In order to compute Cp we 
rewrite its formula in terms of the solutions of the famous Lane-Emden ODE 
with standard initial data. In the usual notation (as in |Ch67j ) this ODE is 

0. 

1, (106) 
0. 

The solution to this ODE for a particular choice of n is often referred to as 
the standard polytrope of index n. It is well-known (c.f. [Ch67] ) that the 



de e ^ n 

0n(O) 
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standard polytropes for n G [0, 5) first cross the £-axis at a finite distance 
from the origin. This first zero is often denoted £ n . 
Explicit solutions are only known for three indices: 

0o (0 = (107) 

em = (108) 

e 5 (0 = jl — , (109) 



giving £ = \/6, £l = 7T) and £ 5 = oo. Of equal importance is the slope of 6 n 
at the first zero. In the cases above, we have: 

* (6) - (no) 

lim^K) = 0. (112) 

We next explore rescaling the standard polytropes in order to find func- 
tions which satisfy our equation for (pp ( ITT]) . We first note that the polytropic 
indices arising in the determination of Cp range over (3, 5] (so that n = 1 is 
clearly avoided in our considerations). We make the following definition (for 
n ± 1) 

7n (0 = a o - 1 ^9 n (A,0 I (H3) 
and note the following consequences: 

(114) 
(115) 
(116) 





= o, 








din (f) 

d£ 


n+l 

— a n A n 






— a n n n 





These change our ODE to 



-In A'i^, O n+l A^, -In 

_ An-l In „ A n-1 l n < „n A n-1 -,n n 1 T\ 

anAn IPi + -TT a nA n — + a n A n 7 n = 0, (117) 
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which reduces to 



d In 2 d'Jn n -x n n /i 1Q\ 

W + n ln = ' ( } 

7n (0) = a-^r 1 , (119) 

din, 



d£ 

We clearly must have 



(0) = 0. (120) 



ntf) = (121) 

where n runs from 5 down to 3 as (3 runs from | up to infinity. Equivalently, 
we have 

77—2 

0(n) = -. (122) 

77 — 3 

It is also clear that we will need 

a nW = c((3)^. (123) 
The determination of A n m) comes from the second boundary condition of 

dln(P) ( £n{{})\ _ x ,^k d6 nm 

(125) 



d£ \A nW J <*>) »W di 



A 2 



e 



n(/9) 



Some algebra reveals that 



««(/?) J 

which leads to the following formula for Rp in terms of the standard Lane- 
Emden data: 

Rfl = (-^)^)(W))) 1-2/3 W))"' 1 - (127) 

This in turn leads to the following very useful (but decidedly more cumber- 
some) formula for Cp: 

_ / / 3 (-^(/ ? )^ (/ 3)(^(/3))) 2/3 " i y 

C "- { (2/3-3)^)(c(/3))^ I • (128) 
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Incidentally, we have the following formula for 



which we do not expand further for reasons of brevity! 

6 Numerical Results 

It will be beneficial to rewrite our formula for Cp in terms of the standard 
polytropic index n((3): 

n(/3)-3 

c = ( w)-2)(-^)^(^))^ V' w " 

We first compute the value of Cp as (3 = 3/2. 
First, recall the following facts ( |Bu78| ): 

lim-eX(^n) = V3, (130) 

n— >b 

lim(5-n)£ n = — . (131) 
Since n(3/2) = 5 and c(3/2) = 87r 2 /15, we see that 

c i = s©'' (132) 

which reproduces the value found in |KTZ08] despite the fact that our anal- 
ysis does not apply to the limiting case of non-compactly supported mini- 
mizers. Note that the polytrope of index 5 (commonly referred to as the 
Plummer Sphere in the astrophysical literature) is not compactly supported. 

Using Maple to run the numerical approximations for the standard poly- 
tropes yields the following plot of Cp (displayed with the bounds found 
in [KTZ08J and a vertical line at /3 = 3/2 indicating that Cp = for 
1 < < 3/2): 
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Incidentally, we can also give an improved upper bound over that listed in 
|KTZ 08j. To see this, note that the upper bound in (jSJ) is convex decreasing 
between j3 = 3/2 and /3 ~ 3.6649 beyond which it is strictly increasing for all 
j3 and converges to a finite value (namely 45/87T 2 ) as tends to infinity. In 
contrast, the numerical evaluation of Cp is a decreasing function of (3. Hence 
we can improve the upper bound given in (jSJ) by simply replacing it with 
its convex hull. We find (using Maple to estimate the minimum) that the 
improved upper bound takes the constant value 0.20269 for /3 > 3.6649. 
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7 Asymptotics 



For the limiting behavior as (3 tends to infinity, we need that 

167T 2 

lim c {p) = — . (133) 
We also note that since Hindoo n(/3) = 3, the only terms that contribute are: 

p _ 3 (~^3^3fe)) 2 /-, C XA\ 

Loo ~ 16vr2 • [L6V 

This exact expression is a little less than illuminating since 6% is not known 
explicitly However, there are extensive numerical data available. Referring 
to [Ho86l p. 407], we see that -£|0£(£ 3 ) « 2.018236. Thus, we can at least 
conclude that 

Coo « 0.077383. (135) 

In comparison, Theorem I of |GS85j (recalling that the total mass is 1 
in our considerations) requires that initial data have £°°-norm less than 
40 -3 ~ 0.00002. Of course, Glassey and Schaeffer did not aim for the op- 
timal constant and so were generous in their estimates. In comparison, the 
lower bound given in |KTZ08j is approximately 0.049438. 

Since the standard polytrope of index 5 is known explicitly, we can find 
an asymptotic expression for Cp when f3 is sufficiently close to 3/2. We begin 
by examining an identity involving the zeroes of the standard polytrope of 
index n: 

n + l Jt {0n{r)r l r 2 dr 

(5-n)6. (-£X(U) 2 ' 1 } 

(this is essentially a reformulation of the identity found in l4.1l for the standard 
polytropes). The right-hand side limits to a finite value as n approaches 5 
([Bu78]): 

fw = ^r- (137) 

So, for n sufficiently close to 5, we have that 

{. „ 16 <" + 1 ) , (138) 
7rv3(5 — n) 
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and accordingly, for j3 sufficiently close to 



16 f 4£-3\ / 327r 2 (/3-l) 3 
p ~ 3tt V2/3 - 3 J V3/3(2/5 - l)(3/5 - 2) 

Finally, this yields an asymptotic expression for C« near |: 



3tt 



/3 



16 \ 4/3-3 



v f 3/3(2/3- 1) (3/3-2) 
V 32vr 2 (/3- l) 3 



(139) 



(140) 



Plot of and Approximation (140) 
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